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Abstract 

In these lectures I examine some of the principal issues in cosmology 
from a particle physics point of view. I begin with nucleosynthesis and 
show how the primordial abundance of the light elements can help fix the 
number of (light) neutrino species and determine the ratio rj of baryons to 
photon in the universe now. The value of g obtained highlights two of the 
big open problems of cosmology: the presence of dark matter and the need 
for baryogenesis. After discussing the distinction between hot and cold 
dark matter, I examine the constraints on, and prospect for, neutrinos as 
hot dark matter candidates. I show next that supersymmetry provides a 
variety of possibilities for dark matter, with neutralinos being excellent 
candidates for cold dark matter and gravitinos, in some scenarios, possi¬ 
bly providing some form of warm dark matter. After discussing axions as 
another cold dark matter candidate, I provide some perspectives on the 
nature of dark matter before turning to baryogenesis. Here I begin by 
outlining the Sakharov conditions for baryogenesis before examining the 
issues and challenges of producing a, large enough, baryon asymmetry at 
the GUT scale. I end my lectures by discussing the Kuzmin-Rubakov- 
Shaposhnikov mechanism and issues associated with electroweak baryo¬ 
genesis. In particular, I emphasize the implications that generating the 
baryon asymmetry at the electroweak scale has for present-day particle 
physics. 


1 Introduction 


There is a symbiotic relationship between particle physics and cosmology. This 
is not surprising since both deal with physics in similar environments. Cosmol¬ 
ogy, the physics of the early Universe, is concerned with matter at the high 
temperatures characterizing the Universe at this epoch. Particle physics, the 
physics of fundamental constituents and their interactions, deals with phenom¬ 
ena at very short distances which can only be probed at high energy. Since 
high temperatures and high energies are synonimous, not surprisingly particle 
physics and cosmology are deeply intertwined. 

If one begins to think of which aspects of particle physics are relevant to 
the early Universe, one arrives soon at a very long list. For convenience, I have 
split up this list into four broad categories. The first category includes what 
might be called Planck scale physics. These are interactions whose natural 
scale is of order of the Planck scalej^ Mp ~ 10^® GeV, and which are of (likely) 
importance in the early Universe. A well known example is provided by Grand 
Unified Theories (GUTs). The second class revolves around the physics of 
light excitations. There are a variety of established, or postulated, nearly 
massless particles like neutrinos and axions, which may well play an impor¬ 
tant role in the Universe’s energy density and could have had a role in creating 
structure. The third category encompasses stable, or long lived, heavy par¬ 
ticles. These particles, of which the LSP of supersymmetric theories is a prime 
candidate, can be part of the constituents that make up the dark matter of 
the Universe. Or, if present in large enough quantities, as in the case of super¬ 
heavy magnetic monopoles, they can have a nefarious role in the evolution of 
the Universe. In the final category, I include the consequences of symmetry 
breakdown. One suspects that phase transitions, of different kinds (e.g. the 
one that gave rise to the inflationary phase of the Universe), play a crucial role 
in the evolution of the Universe. In addition, the breakdown of discrete sym¬ 
metries, like CP, or of some continuous symmetries can influence substantially 
the resulting cosmology. 

Just as different aspects of particle physics affect the evolution of the Uni¬ 
verse, conversely the physics of the early Universe also has an important bearing 
on particle physics. That is, cosmological observations can help inform particle 
theory. For instance, as I will show below, the primordial abundances of light 
elements effectively constrains the number of light neutrino species. Eventually, 
high precision data on the angular and power spectrum of the cosmic microwave 
background radiation should help pin down the neutrino mass spectrum. Sim¬ 
ilarly, as we will see, the precise nature of baryogenesis deeply influences the 

^The Planck mass Mp is the mass scale derived from Newton’s constant. In the natural 

— 1/2 

system of unity we are using, where h = c = k = 1, Mp = ' . 
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view one has of the sources for CP violation. 


In many instances, the back and forth relation between particle physics and 
cosmology has proven very stimulating for both fields. Baryogenesis provides 
perhaps the best example of this symbiotic relationship. The Sakharov condi¬ 
tions for baryogenesis in the Universe, enunciated in 1967[D, were first made 
manifest in GUTs about a decade later and contributed to the enormous interest 
in these theories. However, GUTs also overproduced magnetic monopolesj^ cre¬ 
ating a cosmological crisis which was only resolved through the development of 
the inflationary Universe scenario]^. Although it was pointed out by’t HooftQ 
already in 1976 that, as a result of the chiral anomaly, baryon number is not 
exactly conserved in the Standard Model, the rate for these processes seemed 
insignificantly small to be much more than a curiosity. However, about a decade 
later Kuzmin, Rubakov and ShaposhnikovQ showed that these processes could 
be important at temperatures near the electroweak phase transition, opening up 
the possibility that baryogenesis occurred much later in the evolution of the Uni¬ 
verse than hereto believed. Bounds on the Higgs mass obtained at LEP in the 
1990’s, however, suggested that this interesting cosmological scenario was only 
tenable if there were additional CP violating phases, besides the usual CKM 
phase of the Standard Model. What the next development in this saga will be 
is unclear. Nevertheless, it is obvious that, at least in this area, cosmology and 
particle physics are deeply intertwined. 


2 Primordial Nucleosynthesis and the Number 
of Neutrino Species 


I begin my lectures by discussing nucleosynthesis. Although this material is well 
known|P, its affords me a way to introduce, in a familiar context, a number of 
concepts which will be of use later. Furthermore, nucleosynthesis is also the 
first area where a cosmological observation had a direct bearing on particle 
physics 10, so it makes sense to begin here. 

One has known for a fairly long time that the bulk of the Helium present in 
the Universe is primordial^. Although a small amount of the approximately 
25% mass fraction of Helium was generated in stars, all the rest was generated 
by nucleosynthesis in the early Universe. The calculation of this primordial 
fraction of Helium, Ip, by Wagoner, Fowler, and Hoyle|^ in the late 60’s was 
one of the early triumphs of cosmology and remains an important milestone 
for our understanding of the Universe. When one examines the ingredients 
that lead to a prediction of Ip ^ 0.25, two play a crucial role. These are the 
energy density of the Universe at the time when the neutrons and protons go 
out of equilibrium and the temperature where enough deuterium is formed. The 
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former, in detail depends on the number of light neutrino species N^. The latter 
is related to 77 = nsln-y, the ratio of baryons to photons in the Universe now. As 
we shall see, the ratio rj is an important cosmological parameter, related both 
to the quantity of dark matter in the Universe and to the asymmetry between 
matter and antimatter in the Universe. Nn, on the other hand, is a crucial 
number for particle physics. This quantity is now known to great accuracy as 
a result of precise measurements of the width of the Z boson. However, before 
these measurements N^, already could be determined reasonably well indirectly 
through its numerical influence in predicting the Helium mass fraction kp|Q. 


I will sketch now the calculation of Yp, focusing particularly on how the 
final answer depends on and 77 . The crucial concept to understand is the 
idea of freeze-out, or decoupling, of physical processes in the evolution of the 
UniverseQ. This occurs when the interaction rate T = n{av) for the process 
in question becomes slower than the Universe’s expansion rate. In the standard 
Big Bang cosmology, this latter rate is given by the Hubble parameter H, which 
scales with the Universe’s temperature as 


T2 

Mp 


( 1 ) 


If r for certain processes is much greater than H, then these processes are 
in equilibrium in the Universe. Conversely, if T <C i/, the interaction rate is 
too slow compared to the Universe’s expansion rate to keep these processes in 
equilibrium. The freeze-out temperature is the temperature at which T ~ if. 
That is, it is the temperature (or time) when certain processes begin to go out 
of equilibrium in the Universe. 


There are two important moments for nucleosynthesis. The first of these is 
related to the freeze-out of the weak interactions between neutrons and protons. 
Above this freeze-out temperature, neutrons and protons are in equilibrium 
through the weak interactions 
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The rate for these processes scales as T ~ GpT^, and the ratio of neutrons to 
protons is fixed by their mass difference, through the usual Boltzmann factor 

n/p=e-^™/^. (3) 


Freeze-out occurs when the Universe cools to a temperature of around 10^° °K ~ 
I MeVQ, when T ~ The freeze-out temperature T* fixes the ratio of neutron 
to baryons at that time in terms of the Boltzmann factor: 
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The Helium mass fraction Yp depends on Xn(T*), and the particular value 
of Xn{T*) one obtains depends in detail on This latter assertion is easily 
verified by examining Einstein’s equations for a Friedmann-Robertson Walker 
Universe,^ which relate the Hubble parameter to the matter density: 


= 



( 5 ) 


The matter that drives the expansion of the Universe at this time is composed 
of the states which are still relativistic then: the photons, the electrons and 
positrons, and the N^, species of neutrinos. Thus 

p = Pe±+ Pujv = + ^ + 

Here a is the Stefan-Boltzmann constant and the different weights take into 
account the different statistics between bosons and fermions and the fact that 
neutrinos have only one active helicity component.From (5) and ( 6 ) one sees 
that the expansion rate H at T* scales as i? ~ (22 + 7N^y^'^T*'^/Mp. Since 
T ~ GpT*^, it follows that T* depends on Ni, as 

T* ~ (22 + . (7) 

In view of the above and Eq. (4), one sees that the neutron to baryon ratio at 
freeze-out X„(T*) increases if the number of neutrino species N^, increases. 

After neutron-proton freeze-out, the ratio A„ decreases exponentially be¬ 
cause of neutron decay, so that at any time t after t* 

A„(t) = A„(U)e-^ , (8) 


with Tn being the neutron lifetime. Helium nucleosynthesis occurs at a time 
td (or temperature Td) when enough deuterium is formed {Xd ~ A„), since 
(almost) all deuterium transmutes directly into Helium through the reaction 
d + d ^ He -I- 7 . Because the reaction n + p ^ d + j has a fast rate, the 
deuterium fraction is fixed by a Boltzmann factor. One has 
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In the above B is the deuterium binding energy, B ~ 2.2 MeV, and nB{T) is 
the density of baryons at the temperature T. Nucleosynthesis starts at T = 

^The value given in Eq. (4) corresponds to that obtained for Nu = 3. 

^In Eq. (5), R is the scale parameter characterizing the FRW Universe.The curvature 
term, at this early stage of the Universe, can be safely neglected and is omitted from this 
equation. 
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Figure 1: Predicted values for Yp for different values of Ni, 


when the ratio above is of 0(1). Thus, the temperature Td is intimately related 
to the baryon number density at that stage of the Universe. In turn, one can 
relate nB(Td) to the density of photons at that epoch n-y{Td) which just depends 
on Td, n-y(Td) Tj. The argument is simple. Because the ratio of the baryon 
to photon densities is independent of temperature, knowing the baryon to 
photon ratio now, rj, it follows that 

«B(Td) = WjiTd) 7?T| . (10) 


Numerically, one hndsj^ that Td 10® °K ^ 0.1 MeV. One sees from Eqs. 
(9) and (10) that if the ratio 77 decreases, so does the temperature Td when 
nucleosynthesis starts. Basically, for smaller rj one needs a larger Boltzmann 
factor, At Td, because of neutron decays, the ratio Xn{Td) — 0 . 12 , 

roughly half of what it was at freeze-out. Since at Td, Xn ~ Xd and all the 
deuterium is transmuted into Helium, one expects 

Up ~ 2Xr,iTd) ~ 0.24 . (11) 


Precise results for the Helium mass fraction Ip depend on the actual value 
of Ni, and 77 (as well as on the neutron lifetime, t„). As we saw, XniT*) is larger 
the larger Nn is. Thus Yp increases with increasing Ny. Similarly, a larger 77 
also leads to a larger X„(^) and hence a larger Yp. Fig. 1 shows the results 
of a detailed calculation of the primordial abundance of He, plotted as a 
function of 77 , for three different values of N,, = 2,3,4). Clearly, if one 

knew 77 , from estimates of Yp one could infer a value for One way to infer 
77 is to compute also the primordial abundances of other light elements, besides 
Helium, produced by nucleosynthesis. Demanding concordance of these results 
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fixes 77 and one can then infer a value for Nu from cosmology Q. The results 
of the Chicago group [p^ in the early 1980’s, shown in Fig. 2, using a range 
0.22 < Ip < 0.26 for the Helium abundance, predicted 

iV, < 4 . (12) 


This cosmological inference was verified at LEP and the SLC almost a decade 
later, by studying e+e" scattering at the energy of the Z boson mass. The 
amplitude for the annihilation of and e~ into a fermion-antifermion pair 
depends on the width of the Z 


A{e+e //) 
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(13) 


This width, in turn, depends on the number of light neutrinos species—where 
light here means ttIi, ^ One has: 


Tz = T{Z charged states) -|- NvT{Z vv) . (14) 

Fig. 3 plots some early data from the ALEPH|1^ collaboration at LEP which 
clearly shows that N^, = 3 is preferred. The most recent compilation of results 
from all the four LEP collaborations |T^ doing a Standard Model fit, gives 

= 2.993 ±0.011 . (15) 


A less accurate, but more direct measurement of the, so-called, invisible width 
of the Z —assuming that this width is due to Z decays into neutrino pairs—gives 
instead 

= 3.09 ±0.13 . (16) 

Thus, there is now strong evidence that Ni, = 3. 


Given these results, the present-day discussions of nucleosynthesis take Ni, = 
3 and try to get stronger limits on the baryon to photon ratio rj from the demand 
of concordance of all the primordial abundances. A recent example of such an 
analysis is the work of Copi, Schramm and Turnerj 14 , whose results are depicted 
in Fig. 4, yielding for rj the range 


2.4 X 10"^° < 77 < 4.2 X 10"^° 


(17) 


A recent study of updated data for primordial ^He by Olive, Skillman and 
Steigmanj^ pins down Fp in a narrow range 


Yp = 0.234 ± 0.002 ± 0.005 , 


(18) 
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Figure 4: Concordance of meausured primordial element abundances, from 
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where the first error is statistical and the second is an estimate of the possible 
systematic error. Using the above, the 95% confidence limit for Ip gives Yp < 
0.244, which allows Olive, Skillman and Steigman|^ to set a 95% CL for rj of 

?7 < 3.8 X 10-1° (95% CL) . (19) 

This result is in agreement with the recent work on the primordial abundance of 
deuterium obtained by studying quasi-stellar objects (QSO) absorption linesQ, 
which infers a rather high primordial deuterium abundance. However, very 
recent work by Tytler et aZ.|j^, based on two correlated QSO observations, 
obtains a discordant, very low, primordial deterium abundance yielding large rj 
values: 

5.1 X 10-1° < 77 < 8.2 X 10-1° . (20) 

Such values correspond to a range for Yp (0.246 < Yp < 0.282) above the 95% 
limit of Olive, Skillman and Steigman[p^. 

Clearly, the situation at the moment is still unsettled and it is difficult to 
draw strong inferences. Possibly, the simplest assumption to make is that the 
actual value for Yp is subject to much stronger wstematic uncertainties that 
those assumed by Olive, Skillman and Steigman||^. In what follows, we shall 
take the value obtained by Copi, Schramm and Turner|^ for 77 but, following 
their suggestion, shall boost it to cover a 2cr range. Then one has 

1.9 X 10-1° < ?7 < 5.8 X 10-1° (2cr range) , (21) 

which is a range broad enough to encompass all recent determinations. 


3 Two Open Problems in Cosmology: Dark Mat¬ 
ter and Baryogenesis 


The ratio 77 , which we just saw is important for nucleosynthesis, lies at the heart 
of two of the biggest open problems in cosmology today, those of dark matter 
and of baryogenesis. Recall that 77 was the ratio of the number density of 
baryons to photons in the Universe now. The photon density itself is extremely 
well known from the measurement of the temperature of the cosmic background 
radiation Q 

Ty = (2.726 ± 0.005)°K , (22) 


yielding 


[ q^dq 1 

J TT^ (e'J/'U — 1) 


[0.625 ~ 400 cm-° . 


(23) 
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Therefore a value for rj serves to fix the energy density of baryons in the Universe 
now: 

PB = ruNriB = mNVn^ , (24) 

where iriN is the nucleon mass. Using Eqs. (21) and (23) one finds 

1.3 X 10-3^ g/cin < PB < 4 X 10-31 g/cm^ . (25) 

This value is interesting since, as we shall see below, it is a few percent of what 
is needed to close the Universe. 


If one does not neglect the curvature term, Einstein’s equations in a Fried¬ 
mann Robertson Walker Universe have the form 


/ r\ SttGn k 

\r) ^ 3 ■ 


(26) 


The constant k here describes the geometry of the Universe. If fc > 0, the 
Universe is closed. If fc < 0, it is open. Finally, if fc vanishes, one has a Universe 
with no curvature—a flat Universe. It is useful to consider the quantity n(t), 
which is essentially the ratio of the matter density to the square of the Hubble 
parameter 


H(f) 
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Using Einstein’s equations, one sees that H(t) characterizes the geometry 
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with H > 1 corresponding to a closed Universe, and H < 1 corresponding to 
an open Universe. The value of H(t) at the present time Ho depends on how p 
compares to the, so-called, critical density 


3gg 

SttGn 


(29) 


with Ho being the value of the Hubble parameter now—the Hubble constant. 


From Eq. (28) one sees that a flat Universe has H(f) = 1. Thus the ratio of 
the Universe’s density now to the critical density: 

Ho = -^ (30) 

Po 

directly informs one about the Universe’s geometry, with Hq = 1 corresponding 
to a flat Universe. Unfortunately, the Hubble constant itself is not that easily 
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determined. Conventionally, one writes:^ 


Ho = 100 


h 


Km 

Mpc sec ’ 


(31) 


and one typifies the uncertainty in Ho through a range for h. Traditionally, this 
uncertainty corresponds to h lying in the range 

0.5 < h < 1 , (32) 

although a more modern determination |l^ gives 

/i = 0.6 ±0.1. (33) 


Numerically, one finds that the critical density has the value 

Po = 1.9 X 10“^^ h^ g/cm^ . (34) 

If the density of the Universe is above po the Universe is closed. Clearly, if 
baryons dominate the energy density of the Universe, Eq. (25) tells us that the 
Universe is open. 

If we denote the baryonic contribution to Uq by Ub = Pb/Po, using Eq. (25) 
and (34), one has|l4| 

0.007 < UBh^ < 0.021 . (35) 

This equation is remarkable in several ways. First, it appears that Ub itself 
is much bigger than the value one would infer from the amount of luminous 
matter in the Universe. Using h = 0.6 ± 0.1, from Eq. (35) one sees that Ob 
ranges from about 0.014 to 0.084. On the other hand, the best estimates of the 
fraction of luminous matter in the Universe ]^ give a range 

Oluminous ^ 0.003 - 0.017 , (36) 

half an order of magnitude smaller. So one infers that there is substantial non- 
luminous baryonic dark matter. The existence of this dark matter is also inferred 
from the observed flat rotation curves in spiral galaxies. Normally, outside the 
luminous body of the galaxy, one would expect the circular velocity to drop 
as but it does not, as shown in Fig. 5. From these measurements, one 

deduces values pl|| 

flrot. curves ^ 0.03-0.10 , (37) 

much more comparable to those for Ob. 

Mega parsec (Mpc) is 3 X 10® light years. 
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Second, since fie is much bigger than the contribution to the energy density 
made by photons, neutrinos and electrons,^ if flo — then there is an enor¬ 
mous fine-tuning problem. In this case, the parameter X{t) now, Xo, is roughly 
of order ten to one hundred: 


^ " Ob ■ 


(38) 


However, X{t) being the ratio of the curvature term to the energy density term 
[cf. Eq. (28)] scales as| X{t) ~ R^{t) ~ T-^. Hence X{t) ~ -1/Ob (To/T)^ 
and therefore 

0(r) ~ ^ . (39) 


To get Oo ~ Ob now, in the early Universe the density must have been unbe¬ 
lievably close to the critical density. For instance at the Planck temperature, 
Tp ~ 10^2 °K, 0(Tp) ~ 1 - 0(10-62) ! 


The solution to the fine-tuning problem above is provided by inflation. In 
an inflationary Universe, there is an exponential growth of the scale factor at 
early times. Effectively then the curvature term k/R^ is totally negligible and 
the latter evolution of the Universe corresponds to that of a fiat Universe, with 
keff = 0. Thus, if one wants to avoid fine-tuning as a result of inflation, then 
H(t) = 1 and the Universe is always at the critical density. In this case, the 
value of Hb obtained, since it is in the percent range, tells us that the Universe 
is dominated by non-baryonic dark matter. 


Besides the theoretical bias for considering H = 1, there is actually obser¬ 
vational evidence for H being greater than Hb, obtained by reconstructing the 
energy density from the flow of peculiar velocities in superclusters of galaxies. 
It appears that the values of Ho one infers are largest when one measures the 
density on the largest structures in the Universe, as shown in Fig. 6. All the 
data on Ho has been summarized recently by Dekel, Burnstein and White 
who, if one assumes that there is no cosmological constant, give the following 
range for this quantity: 

0.3 < Ho < 1.3 . (40) 

The lower bound above comes from the cosmic velocity flows, while the upper 
bound comes from the age of the Universe (assuming h = 0.6 ± 0.1). Fig. 7 
summarizes these results, allowing for the possibility of a cosmological constant. 

®The energy density of photons and neutrinos follows directly from the temperature of the 
cosmic microwave background radiation T^, with p-y = aT/ and T^, = due to 

photon reheating. H Charge neutrality requires rie = rip and hence, because of the proton- 
electron mass difference, the energy density of electrons in the Universe is negligible. 

®I neglect here, for simplicity, the fact that matter dominates over radiation in the lat¬ 
ter stages of the Universe. This changes the fine-tuning problem only qualitatively, not 
quantitatively. 
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Figure 6: Variation in SIq as a function of the scale of the structures meausured 


Very recent data on type-I supernovas at high redshift has pro¬ 

vided some evidence, at the 3cr level, that the Universe’s expansion is actually 
accelerating rather than decelerating. Since the deceleration parameter 
meausures the difference between the contribution of matter to U and that of 
the cosmological constant 

= (41) 

this data suggests a non vanishing cosmological constant. However, these results 
are based on the assumption that supernovas are standard candles. Further, 
they are critically dependent on the highest redshift supernovas observed. There 
is also some disagreement among the results of the two groups. For a flat 
Universe, |2^ obtains Um = 0.32 ± 0.1 and Ua = 0.68 ± 0.1, while Q find 
Um = 0.6 ± 0.2 and Ua = 0.4 ± 0.2. In my view, it is probably to early to 
abandon the idea of a Universe where only matter (of all types) contributes to 
give U = 1. However, these results give one pause. 

The parameter ry, besides fixing Ub and adumbrating the dark matter prob¬ 
lem, has another role, ry is also a measure of the amount of matter-antimatter 
asymmetry in the Universe. From observation, it appears that the Universe is 
matter dominated, with little or no antimatter, The observed antiprotons 
in cosmic rays, whose typical ratio to protons is of 0{p/p ~ 10“^), are entirely 
consistent with the flux coming from pair production. Furthermore, no charac¬ 
teristic 7 -rays are seen in the sky which could arise from p — p annihilations. 
If the Universe had islands of antimatter, one would expect such signals to be 
present. In addition, there are also theoretical difficulties in assuming that the 
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Universe was matter-antimatter symmetric in its late evolution. In this case, one 
can estimate the amount of matter that would remain after the p and p in the 
Universe go out of equilibrium around T ~ 0(1 GeV). Below this temperature 
inverse annihilations (27 —> p+p) are blocked and the direct process p-\-p 27 

considerably reduces the number of protons (and antiprotons) compared to that 
of photons to values of 77 ^ 10 “^® |^. 

For these reasons, 77 ^ 0(10“^°), as observed, is evidence that there was 
some primordial baryon asymmetry. That is, really. 


(rtB - »b) 

n-y 


(42) 


If the value of 77 were to codify an initial asymmetry for the Universe, this would 
appear to be a pretty mysterious initial condition. Fortunately, as Sakharov Q 
first pointed out, it is possible to generate such a baryon-antibaryon asymmetry 
dynamically and so 77 can be a reflection of some primordial processes. To 
generate such an asymmetry dynamically, as we will discuss later on in much 
greater detail, the underlying theory must violate baryon number, as well as 
C and CP. Thus it appears that even though 77 is an important cosmological 
parameter, its origins are tied to particle physics also! 


4 Hot and Cold Dark Matter 


An important classification scheme for dark matter is whether the relic dark 
matter candidates were created by a thermal process or as a result of some 
non-thermal process (e.g. in a phase transition). Thermal relics can be further 
distinguished by whether they were relativistic or non-relativistic at the time 
their interaction rate fell below the Universe’s expansion rate. Relics which were 
relativistic at freeze-out are labeled hot dark matter (HDM), while relics which 
were non-relativistic at freeze-out are called cold dark matter (CDM).[[] 

Particle physics provides possible dark matter candidates in all these cate¬ 
gories. Neutrinos, neutralinos and gravitinos are thermal relics, while axions are 
an example of a non-thermal relic. Neutrinos are a prototypical hot dark matter 
relic. Neutralinos are an example of cold dark matter, while gravitinos are warm 
dark matter candidates. Because only zero momentum axions can contribute 
substantially to the Universe’s energy density, axions are also cold dark matter 
candidates. In what follows, I will describe some of the characteristics of these 
possible dark matter candidates. 

^There is also warm dark matter (WDM). These are relics which, while relativistic at 
freeze-out, have much weaker interaction rates and so, in some sense, have also some of the 
characteristics of CDM. 
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The contribution of thermal relics to depends on what their abundance was 
when their interaction rate fell below the Universe’s expansion rate. Freeze-out 
occurs when, for the relic x, Tx — For hot relics the freeze-out temperature 
is much greater than the mass of the relic: ^ m^. In this case, at freeze-out 

^ n-y. Because the density of the relic to that of photons is an invariant, the 
contribution to flo of Einy hot dark matter relic is just a function of its mass 
(and T° ~ 3°K). Calling this contribution one finds 


rixpDM] = 


rriyn 


X"‘x 


m^n. 




1 


rriy 


* 192 eVJ 


1 


(43) 


po 

Here g* counts the effective degrees of freedom at freeze-out, with g* = 1 for 
neutrinos. The above shows that particles with eV masses can be cosmologically 
significant, contributing substantially to the Universe’s energy density. This 
observation was first made about 25 years ago by Cowsik and McClelland 
and Marx and Szalay|M with regards to neutrinos. 


Cold relics, on the other hand, undergo freeze-out at temperatures much 
less than their mass: In this case their density is suppressed relative 

to the photon density by a Boltzmann factor, so that at freeze-out <C Uy. 
This density, however, can be deduced from the freeze-out condition itself 


Ty^ = nx{(Jv)x 


= H^ l.l{g*Y/^Tl/Mp , 

freeze-out 


(44) 


where g* is the effective number of degrees of freedom at freeze-out and (crw)^ 
is the, thermally-averaged, annihilation rate for the cold dark matter relic x- In 
this case, the contribution to Hq of the relic depends both on this annihilation 
rate and on the ratio One finds, approximately, 1^ 
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H,[CDM] ^ 

a formula first deduced by Zeldovichp^. For a typical ratio rriy^/Ty, ^ 20 one 
needs cross sections of 0{a ^ 10“^® cim). These cross sections are of the typical 
strength of weak interaction processes. It is clearly intriguing that such cross 
sections could have cosmological significance! 


10 ^'cm®/sec 


XI freeze-out 


(45) 


There are no simple formulas to describe the contributions to H of non- 
thermal relics, since these contributions depend in detail on the dynamics. In 
general, for non-thermal relics, the interactions are so feeble that one has always 
T^ iJ; that is, the relics are never in thermal equilibrium. For example, 
as we shall see later on, axions with ma — 10“® eV can close the Universe 
(Haxions — !)■ This means that the number density of axions in this case is 
about 10^ times what it would be if axions were thermal relics (i.e. had a 3° K 
temperature). So, if axions are the dark matter in the Universe, they obviously 
had a highly non-thermal origin. 
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5 Prospects of Neutrinos as Dark Matter 


Neutrinos are interesting candidates for dark matter since their properties fit the 
required profile. Furthermore, neutrinos are the only dark matter candidates 
whose existence is confirmed experimentally! Originally, neutrinos were thought 
to provide possible examples for both hot dark matter and cold dark matter. 
Because of LEP, we know now that they can only be HDM candidates. Let me 
elaborate on this point. 


Experimentally, one has evidence that the three known neutrinos, and 

Vt, are quite light, with direct bounds on their masses given by[^ 

< 15 eV ; < 170 keV ; < 24 MeV . (46) 


Because the freeze-out temperature for neutrinos is of order of Ty ~ 1 MeV, so 
as not to overdose the Universe and Vr must have masses much below these 
bounds. Hence, it is perfectly conceivable that the known neutrinos are the hot 
dark matter, contributing to Ho an amount 




1 


92 eV 


■ 


(47) 


If heavy neutrinos existed, with masses ^ Tf ^ \ MeV, they could 
be cold dark matter candidates because they have weak interactions. It is 
straightforward, knowing the interaction rate of neutrinos, to calculate their 
contribution to H as a function of the neutrino mass. The result is displayed in 
Fig. 8, taken from^. This figure shows that grows with rrii, up to around 
the freeze-out temperature Tf and then decreases rather rapidly for neutrino 
masses beyond this temperature. However, the existence of further neutrino 
species with rrii, < \Mz is now excluded by measurements of the Z-width at 
LEP, which, as we saw earlier, gives = 3 to high accuracy. Thus the window 
for heavy neutrino CDM is closed. 

Although neutrinos are not cold dark matter candidates, they remain excel¬ 
lent prospects for hot dark matter. Nevertheless, because one knows that hot 
dark matter alone cannot describe the power spectrum of density fluctuations!^ F| 
one expects Hy < 1. Acceptable fits to the power spectrum of density fluctua¬ 
tions suggest typically Hy ~ 0.2. Using /i^ ~ 0.3, if neutrinos are the HDM, 
this gives 

~ 5 — 6 eV . (48) 

This ratio is consistent with the bounds p0| given in Eq. (46). Unfortunately, 
however, there is as yet no direct particle physics evidence that the known 

®HDM neutrinos, because they are so light, have a free streaming length of 0(10 Mpc). 
As a result, they cannot account for the formation of structure at small scales. 
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Figure 8: Contribution of neutrinos of different masses to from 


§ 


neutrinos have masses that satisfy Eq. (48). Nevertheless, there is tantalizing 
indirect evidence for neutrino masses (through hints of neutrino oscillations) 
and this evidence is compatible with Eq. (48). Because of its importance to the 
issue at hand, I review next some of this information and its implications. 

First, let me make a comment on prospects for improving the direct neutrino 
mass bounds quoted in[pO[. Clearly, even though the kinematical techniques that 
give the bounds on and can perhaps be improved somewhat, there is 
no hope to directly measure masses in the few eV range for these particles. 
This is not so for r'e- In fact, the tritium /3-decay experiments that lead to the 
bound of < 15 eV quoted by the Particle Data Group||3^, actually all have 
sensitivities of order 1-2 eV! The reason for the much weaker bound quoted, 
is that all the latest precision experiments are plagued by an anomalous 
unexplained excess of events beyond the tritium beta decay endpoint. This 
excess actually leads to an average mass-squared that is negative {{mlj = 
(—27 ± 20) eV^).Q Until this excess is understood, one cannot set a real 
bound for , although the potential sensitivity to eV masses is there. 

In this context, one should mention a different piece of evidence that suggests 
that i^e itself cannot be the dominant form of the HDM. This latter constraint 
comes from double-beta decay, where searches for the neutrinoless mode in ^®Ge 
decay^^ lead to a limit on the effective neutrino Majorana mass responsible 
for this process of [| 

{m^)ee < (0.5 - 1.5) eV . (49) 

®The uncertainty in the bound of Eq. (49) reflects an uncertainty in the calculation of the 
relevant nuclear matrix element. 
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Here {mi,)ee is the sum of the neutrino masses entering in the process, convo¬ 
luted with the appropriate mixing matrix element coupling these neutrinos (and 
antineutrinos) to electrons: 


{mv)ee = ^ UliTUi . (50) 

i 

If mixing of electrons to neutrinos other than is not large {Uei — 1), then 
Eq. (49) is also a bound on However, this is not a true bound because if 

neutrinos are Dirac particles, the particle and antiparticle contributions in Eq. 
(50) automatically cancel and (m^)ee = 0. 

Eortunately, one can probe neutrino masses indirectly by looking for evi¬ 
dence for oscillations of one neutrino species into another. If neutrinos have 
mass, neutrinos can mix with one another and this mixing can be revealed 
through neutrino oscillation experiments. At present there are a number of tan¬ 
talizing hints arising from experiments looking for neutrino oscillations which 
have an important bearing on the question of neutrino mass. To discuss these 
experiments, it is necessary first to briefly discuss a bit of phenomenology. 

If neutrinos have mass, the weak interaction eigenstates (the neutrinos pro¬ 
duced by weak interaction processes-e.g. e+z/g), are not the same as the 

mass eigenstates (i.e. the observed particles of well defined mass, denoted here 
by Vi). However, these states are related by a unitary transformation, so that 
each VI {(. = e, r} can be written as a superposition of the Vi\ 

Vi = '^ UiiVi . (51) 

i 

Conventionally, one only examines the 2x2 case, assuming that, as in the quark 
case, the 3x3 matrix will be nearly diagonal with dominant mixing among pairs 
of neutrinos. In this case, for instance, Eq. (51) for v^ and just involves a 
simple orthogonal 2x2 matrix: 

f >ie \ _ f cos9 sin6> \ f \ 

\ Vfj, J ~ y -sinO cos9 J \ V 2 J ' 

The mass eigenstates Vi have the usual quantum mechanical evolution with time: 

kz(t)) = e-*^->,(0)) . (53) 

Imagine then producing at t = 0 a from a weak decay 

|z^(0)) = \ve) = cos0|z/i(O)) -I- sin0|z/2(O)) . (54) 

I am being a little sloppy here not distinguishing between weak interaction eigenstates 
and mass eigenstates. Also, I am explicitly assuming that neutrinos are Majorana particles. 
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At a later time, because the states vi and V 2 have different masses, this state will 
evolve into a superposition of both \ve) and \v^) = — sin0|j^i(O)) + cos0|z^2(O)). 
That is 


\v{t)) = cosd|:/i(t)) + sin0|i/2(i)) 

= cos6»e"*^i>i(0))+sin6»e"*^^>2(0)) . (55) 


Using the above, it is easy to calculate the transition probability that an 
initial Ve state has oscillated into a v^i state after a time t: 

= i sin^ 2d[l — cos(Ui — E 2 )t] . (56) 

2 

In all cases of interest |p| :g> rrii. Hence Ei ~ \p\ + with |p| = E^ being 
essentially the neutrino energy. Also, in this case, the time t in Eq. (56) can just 
be replaced by the distance travelled (in units of the speed of light): t = L/[c\. 
Whence, one finds the following formula for the probability that, as a result of 
neutrino mixing, an initial Ve of energy Ei, has oscillated after a distance L into 
a v^-. 


P{ve^Vfj,;L) = sin^20sin^ 


= sin^ 29 sin^ 


{m\ — to|)L 


1.27 


AE^ 
Am^{eV)L{m) 


E^(MeV) 


Of course, 


Pii 


p,Z/) — 1 P{u^ > . 


(57) 

(58) 


One sees from the above formulas that the probability of oscillation is sen¬ 
sitive to the mixing angle 9. Further, if one wants to probe a particular Am^ 
range, then for a given neutrino energy E^, there are appropriate distances L 
where the effect is maximum. If the neutrinos are not nearly degenerate, then 
the Am^ range one wants to probe to find out whether neutrinos contribute sig¬ 
nificantly to the dark matter problem (cosmologically significant neutrinos) is 
Arn? ~ 25 eV^. This is the goal of the CHORUS and NOMAD experiments at 
CERN, which for Am? in this range hope to be sensitive to ^ oscillations 
as low as sin^ 20^,- > 10“^. 

Up to now the CERN experiments have only given limits. |Q However, in 
other regions of Am^ there are various hints of neutrino oscillations. In fact, 
there is an embarrassment of riches! The LSND experiment)^ sees a signal 
of > De oscillations in a narrow range in Am^ around (0.2-2) eV^, with 
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sin 20 ~ 10“^ eV . The ratio of to Ve neutrinos, observed in large un¬ 
derground experiments, arising from decay processes in the atmosphere shows 
a deficit 1^ from what is expected. This atmospheric anomaly can be inter¬ 
preted either as being due to ^ Vt or ^ z/g oscillations, with Am^ ~ 
(0.3 — 3) X 10“^ eV^ and large mixing angles sin^ 20 ~ 0(1).^ Finally, exper¬ 
iments measuring the flux of solar neutrinos, also show a dearth of neutrinos 
compared to the predictions of the, so called, standard solar model. One can 
reconcile the observations of all of these solar neutrino experiments by appealing 
to Vf, —> or Vf. —!■ Vt neutrino oscillations, which are enhanced in matter by 

the so-called MSW mechanism,!^ provided that Am^ ... I'n ■?_ i ow in-5 , 
with rather small mixing: sin^ 20 ~ (4 — 10) x 


(0.3-1.2) X 10-5 eV 


Because we know of only three neutrino species, these hints cannot all be 
true, since we have at most only two mass differences Even if one were to 
eliminate one of the hints (LSND perhaps-since, as Fig. 9 make clear, the 
allowed region is almost ruled out by other negative findings), because the fa¬ 
vored mass differences are small, it appears that to have cosmologically signifi¬ 
cant neutrinos one must have near mass degeneracy. For example, the pattern 
TOj/j ~ ~ ~ 1 —2 eV, with ~ eV^; ^ IQ-^ eV^ would 

explain the solar and atmospheric anomaly. If this were really the case, then 
cosmologically significant neutrinos would produce no signal in CHORUS and 
NOMAD !0 Hopefully, in the next hve years with upcoming neutrino oscilla¬ 
tion experiments (as well, perhaps, with some clarification in the direct tritium 
/3-decay experiments) one should be able to sort out this somewhat confusing 
situation, thereby arriving at a better understanding of whether or not neutrinos 
can contribute to the dark matter in the Universe.]^ 


Before closing this section, I would like to make an important theoretical 
point. Neutrino masses in the eV and sub-eV ranges are very interesting from 
a particle physics point of view, since they are most likely a signal for a new 
large mass scale. In general, because neutrinos are neutral, they can have both 
Dirac - particle/anti-particle—and Majorana - particle/particle masses. That 
is, one can write 

-Cinass = C'z/L + , (59) 

^^Very recent data from SuperKamiokande |^] f avo rs a lower range for Am.^ ~ (0.1 — 1) x 
10~^, while the negative results from the Chooz ]38|| reactor experiment now excludes the 
Ve. oscillation option for explaining the atmospheric anomaly. 

^^It is possible not to discard any experimental hints if one assumes that, in addition to 
i/e, and Vt-, there is an extra sterile neutrino Vs- Then one of the experimental results-the 
solar anomaly-can be interpreted as a ^ oscillation[|l2|. 

^^In this scenario, one has to worry about the double-beta decay limits, since these provide 
effective electron neutrino masses (mu)ee precisely in this range. 

^“^The study of the angular power spectrum of the cosmic background radiation can also 
provide information on this issue, as massive neutrinos can affect this spectrum differently 
depending on their mass.||43| 
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Figure 9: Signal region for the LSND experiment, along with exclusion limits 
from other experiments 
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where C is a charge conjugation matrix. Q Because vl is part of an SU(2) 
doublet, while vb. (if it exists!) is part of an SU(2) singlet, it is clear that 
in the standard model m^,mD and TOmj respectively, carry effective SU(2) 
quantum numbers of 0, 1/2 and 1. In particular, while can be a totally 
independent mass parameter, m d and m\^ must be proportional to the vacuum 
expectation value of an SU(2) doublet and triplet held, respectively. We know, 
as a result of the experimentally very successful interrelation between and 
Mz'- cos^ 0VV = that what causes the breakdown of the electroweak 

theory through its VEV transforms dominantly as an SU(2) doublet. Thus, if 
vb exists, we expect rriD, with rriD ^ m (—the mass of the corresponding 

lepton. 

It was realized long ago by Yanagida[^ and Gell-Mann, Ramond and Slansky, 
that if the Majorana mass of the right-handed neutrinos is very large, 
^ niD, the above scenario produces very tiny neutrino masses. If one 
neglects altogether ~ 0, one has a 2 x 2 neutrino mass matrix of the 

form 

M = ( ° "‘g ) . (60) 

mo rn‘^ J 

This matrix has a very heavy neutrino, mostly vb, with mass and a very 
light neutrino, mostly with mass 


rriy 


ml) 




(61) 


This, so-called, see-saw mechanism Q can produce eV neutrino masses pro¬ 
vided is sufficiently large (e.g. for one has eV neutrinos associated 
with the tau if ^ 10^° GeV!). Thus detecting light neutrino masses is tan¬ 
tamount to discovering a large scale—the scale responsible for the right-handed 
neutrino Majorana mass.|^ 


6 Supersymmetric Candidates for Dark Matter 


Supersymmetric extensions of the Standard Model for the strong and elec¬ 
troweak interactions provide excellent candidates for cold dark matter. Super- 
symmetry, as is well knownis a fermion-boson symmetry. Thus, if it were 
a true symmetry of nature, we would expect a doubling of all the degrees of 

should comment that even if there were no z/r— something I consider unlikely—the 
presence of eV neutrino masses again, most likely, reflects another large mass scale. For 
instance, without uji one can get a Majorana mass for t)y using a doublet Higgs field twice to 
make a triplet. Sucb_interactions are non-renormalizable, but could arise effectively from some 
GUT interactions pTj and are scaled by l/AfQUT* formula like ($)^/AfGUT 5 

with (^) ~ 300 GeV gives eV neutrino masses for Mqut 10^^ GeV. 
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freedom.|^ Thus, in these extensions of the Standard Model, there is a plethora 
of undiscovered particles. Some of these particles turn out to be good dark 
matter candidates. This is plausible because a supersymmetric extension of the 
standard model preserves the strength of the couplings. For instance, the su¬ 
persymmetric vertex joining a squark (the scalar partner of the quark), a quark 
and a gaugino (the spin-1/2 partner of a gauge boson) has the same strength 
as the quark-quark-gauge boson vertex. As a result, (some) of the interaction 
cross sections for supersymmetric (SUSY) particles will have the strength of the 
weak interactions (provided these particles are not much heavier than the weak 
bosons) and thus will satisfy Zeldovich’s criteria for cold relics, Eq. (45). 

Most supersymmetric extensions considered contain a discrete symmetry, 
called i?-parity: 

i? = (-l)3i^+i^+2'S , (62) 

which is -1-1 for particles and -1 for sparticles. If R parity is conserved, then the 
lightest supersymmetric particle, the LSP, by necessity is stable. If it is neutral, 
as is usually assumed to be the case to avoid cosmological difficulties associated 
with their luminosity, the LSP provides an excellent candidate for cold dark 
matter. 

We know that if supersymmetry exists it must be broken in nature. Other¬ 
wise, the masses of the supersymmetric partners, rh, would be the same as that 
of the ordinary particles, m, in gross contradiction with experiment. However, 
we do not know really how supersymmetry breaks down. As a result, which par¬ 
ticle is the LSP is model dependent. Nevertheless, one can make some general 
observations. 

There are three important scales associated with SUSY breaking. The first of 
these is, obviously, the masses of the sparticles themselves, m. The second is the 
scale. A, which is associated with the spontaneous breaking of supersymmetry. 
This is assumed to occur in a, so called, hidden sector, separated from the 
ordinary interactions of particles.^ The last scale, M, is the scale associated 
with whatever phenomena acts as the messenger connecting the hidden sector 
with ordinary matter. This connection is shown pictorially in Fig. 10. Both A 
and M are model dependent, but one expects the masses of the sparticles to be 
of order 

A2 

(63) 

It is an attractive possibility that supersymmetry resolves the naturalness 
problem of electroweak symmetry breaking, related to why the scale of elec- 

^®For technical reasons, connected to the cancellation of anomalies, one needs also to double 
the number of Higgs doublets, as well as provide appropriate fermionic partners to these states. 

^^Separating the process of supersymmetry breaking from ordinary matter is necessary to 
avoid contaminating ordinary matter with interactions we have not yet seen. 
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Hidden Manifest 

Sector Sector 


Figure 10: Scales associated with SUSY breaking 


troweak breaking v ~ 250 GeV is so much less than the Planck mass. For this 
to be the case, SUSY states cannot themselves have masses much bigger than 
V. Hence, it is generally assumed that supersymmetric partners must them¬ 
selves be of mass rh ~ Thus A and M, the parameters associated with 

supersymmetry breaking, are constrained physically to produce h?/M ~ v. 

With this constraint in mind, two main scenarios have emerged, with each 
scenario producing a different LSP. The first scenario arises out of supergravity 
models,!^ where the hidden sector is connected to the ordinary sector by grav¬ 
itational interactions. Here M ~ Mp and thus A ~ 10^^ GeV. In the second 
scenario!^ the messenger sector is associated with gauge interactions with a 
scale around M ~ 10® TeV. To get Is? jM ~ v necessitates then a much lower 
scale A of spontaneous supersymmetry breaking, A ^ 10® TeV. 

In both scenarios the gravitino, the spin-3/2 supersymmetric partner of the 
graviton, becomes massive as a result of the spontaneous breaking of supersym¬ 
metry, with a mass of order 

A2 

■ ( 64 ) 

From the above, one sees that in the supergravity scenario, the gravitino mass 
is also of 0{v )—typical of all the other masses of the supersymmetric partners. 
Thus, in this case, it is generally assumed that the gravitino is not the LSP, 
but that the LSP is a neutralino. This is the lightest spin-1/2 partner of the 
neutral bosonic particles in the theory—the two gauge bosons, 7 and Z, and 
the two neutral Higgs bosons and Hd- In general, the neutralino is some 
particular superposition of all these spin- 1/2 partners 

X = 0^7 -I- azZ + OuiLu + OdiLd , (65) 

where the Oi are model dependent coefficients. In contrast, in the scenario where 
the messenger are gauge interactions, with M ~ 10® TeV and A ~ 10® TeV, the 
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gravitino is extraordinarily light, 


"i 3/2 ^ ^ O(KeV) < m , (66) 

and is the LSP. While neutralinos, with ~ 0(10— 10^ GeV), are typical cold 
dark matter relics, gravitinos, with m ^/2 ~ O(KeV), act as warm dark matter. 

1 discuss both of these cases, in turn. 

The typical supergravity model|^ which gives neutralino CDM is charac¬ 
terized by a set of universal soft supersymmetry breaking parameters, specified 
at the scale A. In addition, it contains as a parameter the vacuum expecta¬ 
tion ratio between the Hu and Higgs bosons: tan/3 = (Hu)/{Hg). The, so 
called, minimal supersymmetric standard model (MSSM)[Q has actually only 

2 additional parameters, besides tan/3, which determine the LSP mass, 
and the coefficients ai in Eq. (65). These are the common mass, TOi/ 2 , of all 
the gauginos and a mass parameter /r characterizing the supersymmetric cou¬ 
pling between the Hu and H^ supermultiplets. However, even in this minimal 
model, the actual contribution of the neutralino LSP to Hcdm depends on the 
neutralino annihilation cross sections. These, in turn, depend on other model 
parameters, the universal soft breaking mass TOq given to all the scalars and 
certain coefficients {A and B) which typify the strength of trilinear and bilinear 
soft interaction terms. 

As a result, even in the MSSM, there is a large region of parameter space 
which produces a neutralino LSP which potentially could be the cold dark mat¬ 
ter in the universe. Typically, what one requires for a viable model is that 
= 0.2 ± 0.1. As can be seen in Fig. 11, there are plenty of models (each 
represented by a dot) which have 50 GeV < < 200 GeV and lead to 

in the desired range, provided that tan /3 is small and the resulting pseudoscalar 
Higgs mass is large {mA 500 GeV).[Q In general, an LSP much 
below about 50 GeV runs into trouble with the negative results from LEP on 
Higgs searches, as well as on the direct production of supersymmetric pairs, 
Thus there are regions in parameter space that are already excluded, serving to 
rule out some potential CDM models. Glearly the discovery of an LSP would 
have an enormous impact on the CDM question, much reducing the parameter 
freedom one still has now, even for the simplest models. 


In gauge mediated supersymmetry breaking models, in contrast, the grav¬ 
itino is the LSP. Here one has much less freedom since there are not that many 
parameters to vary. Gravitino interactions scale as 1/A^, and so are typically 
very much weaker than weak interactions 


0 ' 3/2 


'1 TeV]^ 
A 


^weak ^weak • 


(67) 


28 





CM 



a 



10 10 10 10 

[GeV] 


Figure 11: Contribution to of various neutralino LSP models, from 
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Figure 12: Transfer function for gravitinos and CDM showing their similarity, 
from II 


As a result, the freeze-out of these interactions occurs at an earlier epoch in the 
Universe, when there were more thermal degrees of freedom. Thus, gravitinos 
have a smaller abundance compared to neutrinos of the same mass-H Typically, 
one finds 


03/2/1^ 


r " 13/2 1 

100 

Ll KevJ 

[9*{Tf)\ 


( 68 ) 


where g*{Tf) is the number of degrees of freedom at freeze-out. Because it 
takes gravitinos of mass of order 1 KeV to close the Universe, cosmologically 
significant gravitinos have a smaller Jean’s mass than neutrinos 




M3 


10^^ Mr 


m: 


'3/2 


(69) 


Thus, in contrast to neutrinos, the gravitino free-streaming length is rather 
small, of order A 3/2 ~ 1 Mpc, much closer to that of cold dark matter. Hence, 
gravitinos are typical warm dark matter—matter that is relativistic at decou¬ 
pling but does not form only large structures. Indeed, as I just mentioned, the 
spectrum of density fluctuations for gravitino dark matter is quite similar to 
that of cold dark matter. This is seen clearly in Fig. 12, from a recent study of 
Borgani, Masiero, and Yamaguchi.|| 

Gravitinos, in my view, are not a particularly attractive form of dark matter, 
as to get the needed ^ 3/2 one needs to have the gravitino mass { 1713/2 ~ A^/Mp) 
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Figure 13: Expected rates from neutralino dark matter, from [Q 


finely tuned around a KeV. But this is not the worse trouble! Because of its 
extremely tiny interaction cross section [cf Eq. (67)] gravitino dark matter does 
not have any hope to be detected ever.[^ In contrast, if the CDM is due to a 
neutralino LSP, in principle, it may be detectable by experimental means. 

Calculation of the rates expected in low background experiments (for in¬ 
stance, those using a ^^Ge detector of sufficient mass), depend both on the 
density of LSPs in our galaxy and on the neutralino-nucleon scattering cross 
section. This latter cross section depends again on the various parameters in 
the supersymmetric model. Except for very light nuclei, it turns out that scalar 
exchange dominates, since it leads to coherent scattering of the neutralinos on 
the target nuclei, so that a^A ^ Fig. 13 shows that the expected rates of 
neutralino CDM for a ^^Ge detector are of the order of 10“^ — 10“^ events/Kg- 
day. Given that present-day detectors (e.g. CDMS |^) are operating with at 
best one Kg of Ge, one is still looking for a factor of 10^ — 10^ improvement to 
have some hope of detecting a potential signal for neutralino cold dark matter. 
This is a daunting, but perhaps not impossible, task. As I said earlier, the 
experimental observation of a neutralino LSP in a particle physics experiment 
would give enormous impetus to the lofty goal of direct dark matter detection! 


^®Besides the tiny cross section, KeV gravitinos also give too little energy of recoil. 
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7 Axions as CDM Candidates 


Axions are pseudo Goldstone bosons associated with a spontaneously broken 
global chiral symmetry, [/(1)pq, introduced to “solve” the, so called, strong 
CP problem. The Lagrangian of the electroweak and strong interactions, in 
general, possesses an effective interaction involving the gluon field strengths, 
and their duals G^"^: 


Geff = d^G^GaM.. • (70) 

OTT 

This interaction breaks P, T, and CP and produces a very large neutron elec¬ 
tric dipole moment unless the parameter 9 is very small.^ This is the strong 
CP problem—why is 9 so small? The imposition of an additional global chiral 
symmetry on the standard model suggested by Quinn and myself, [Q essentially 
serves to replace the 9 parameter by a dynamical field—the axion field. In¬ 
stead of the CP violating interaction (70) one now has instead, a CP-conserving 
effective interaction of the axion field a{x) with the gluons: 


>CpQ = y^GrG.^. , (71) 

where / is a scale associated with the spontaneous breakdown of 17(1)pq. 

The axion is the Nambu-Coldstone boson associated with the spontaneous 
breakdown of the t/(l)pQ symmetry. However, because this symmetry has a 
chiral anomaly—reflected in the appearance of the interaction (71)—the axion 
is not truly massless but acquires a small mass.^ This mass is slightly model- 
dependent, but is of order 


ma ^ 


~T~ 


6 X 10®\ 

7 \^]) 


eV . 


(72) 


One sees that, for large /, axions are very light. Since all couplings of the 
axion scales as 1 //, these particles, if they exist, are also very weakly coupled. 
Although axions are not stable since they can decay into two photons, the 
lifetime for the process a ^ 27 scales as t ~ and becomes enormous for 

large /. 


Quinn and I ||^ made the natural assumption that the scale of 17(1)pq break¬ 
ing coincided with the electroweak scale, / ~ u. Unfortunately, these weak-scale 

^®One finds dn ~ ecm]61 _q,nd hence one needs to have 9 < to respect the 

strong experimental bounds on drT] ^ 

^^The interaction (71) produces for the axion field an effective potential which dynamically 
adjusts so as to cancel the 0 parameter. This potential also has a non-vanishing second 
derivative at its minimum, [pl| corresponding to the axion mass. 
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axions have been ruled out experimentally. ||^ If / is not of 0{v), it turns out 
that astrophysics constrains f ^ v. This is easy to understand. Axions provide 
an extremely efficient way to cool down stars, completely affecting their evolu¬ 
tion. Only if / > 5 X 10® GeV, are axion couplings sufficiently weak so as not 
to run into trouble with a variety of astrophysical observations—ranging from 
the evolution of red giants, to properties of the observed neutrino pulses from 
SN 1987a. 1^ For / > 5 x 10® GeV, the axions are so light, so weakly coupled, 
and so long-lived to be effectively “invisible”. |^ However, these invisible axions 
have potential cosmological consequence, and they prove to be interesting cold 
dark matter candidates! Let me review the arguments for this.]^ 

Axions are typical non-thermal relics, since their properties change as the 
Universe evolves. At the [/(1)pq phase transition, which occurs when the 
Universe’s temperature T ~ /, axions are produced as real Nambu-Goldstone 
bosons (rria = 0). At such high temperatures the axion potential due to QGD 
is ineffective and the 9 interaction of Eq. (70) is not cancelled out. As the Uni¬ 
verse cools towards temperatures of order of the QCD-scale Aqqd, T ~ Aqcd, 
two things happen: the axion potential turns on, serving to cancel 9, and the 
axion acquires its mass, which is of 0(Aq(^q//). This relaxation of the axion 
field to its present configuration, however, happens in an oscillatory way. The 
energy density associated with these oscillations, as we shall see, acts as cold 
dark matter, 

The 9 parameter in Eq. (70) can be thought of as an effective VEV for 
the axion field: (a) = 0/, with the correct vacuum state driving (a) —> 0. In 
the early Universe at T ~ /, in this language, the axion field has an effective 
vacuum expectation (a) = 9f. As the temperature lowers towards T ~ Aqcd, 
the QGD potential for the axion turns on and (a) is driven to zero. One can 
study the time evolution of (a) by studying the equation of motion for the axion 
field in the expanding Universe:]^ 


d'^ja) j7(t) d{a) 

df^ R{t) dt 


+ ml{t){a) = 0 . 


(73) 


It is clear from the above that the effect of the expansion of the Universe is to 
provide a drag term for (o). At early times, or high temperatures, the axion 
mass vanishes and (a) is fixed to its initial value (a) = f9. When the axion mass 
ma(t) begins to turn on, as the Universe’s temperature cools towards T ~ Aqcd, 
(a) undergoes damped oscillations about (a) = 0. 


I shall not try to sketch here the computation of the effective energy density 
associated with these oscillations of (a), but refer to Ref.|^ for an elementary 
discussion. I quote, however, the result of a recent detailed calculation @ which 
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gives the contribution to flo of these oscillations. One finds: 


riah^ = C 


f 


1012 GeV 


(74) 


Here C is a constant of 0(1) which depends on the details of the QCD phase 
transition, while the exponent n is near unity, n = 1.18. One sees that if the 
initial value for (a)// = 6 is of 0(1)—as one may expect naively—then these 
oscillations of the axion VEV can close the Universe if / ^ IQi^ GeV. Because 
what is oscillating is (a), these oscillations correspond physically to coherent, 
zero momentum, oscillations of the axion field. Since Pa = 0, axion oscillations 
are prototypical cold dark matter. 

From the above, it appears that coherent axion oscillations can give rise to 
H = 1 provided / ~ 10 ^^ GeV or Wa ~ 6 x 10“® eV. This is predicated on 
having an initial misallignment angle 0 ~ 0(1). However, Linde|^ has argued 
that in inflationary cosmology, with the reheating temperature Treheating < / so 
that there is not a post-inflationary 17(1)pq phase transition, there is no reason 
why the misallignment angle cannot be very small: P 1. In this case one 
could have Ha ~ 0(1) for smaller axion masses (or / ^ 10 ^^ GeV): 

Ha ~ 0(1) if ma ^ 6 X 10-®d^ eV . (75) 


There are other arguments, however, which suggest that axion masses much 
heavier than rria ~ 6 x 10“® eV can close the Universe. These arguments apply 
in inflationary scenarios where the reheating temperature Treheating > /• In 
this case, one must worry about axionic strings formed at the C/(1)pq phase 
transition. The decay of these strings into axions also contributes to the Uni¬ 
verse’s energy density and this contribution can dominate that due to coherent 
axion oscillations. [Q Unfortunately, there is considerable controversy on this 
point, with some authors—notably P. Sikivie and collaborators!^— obtaining 

Ilstring decay ~ ^oscillation; with Others [|7Cf| deducing Hg^rjag decay ^ Hoscillation- 

If one were to believe this latter estimate, then one obtains Ha ~ 0(1) for ax¬ 
ion masses as heavy as rria ~ 10“'^ eV. These masses are perilously close to 
the mass range excluded by astrophysics, rua > 10“^ eV, corresponding to 
/ < 5 X 10^ GeV. 

These controversies may be resolved experimentally if axions are the dark 
matter in the Universe (and hence are also the dominant form of the dark matter 
in our galaxy!). The basic idea for these experiments is due to Sikivie |7l|] and 
uses the fact that axions couple to the electromagnetic field in a way analogous 
to how they couple to gluons [cf. Eq. (71)]: 

-^077 ^ (^^) 
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Figure 14: Results of invisible axion searches 


with 0077 !//• Because of Eq. (76) axions in our galactic halo in the pres¬ 

ence of a strong magnetic field can be resonantly converted into photons in an 
appropriate cavity. Experiments are presently underway at both the Lawrence 
Livermore Laboratory|^ and at Kyoto University!^ which are sensitive to 
“standard” invisible axions if they are the dominant form of dark matter. 0 Fig. 
14 shows recent results from the Livermore experiment [fz^ in the — ma 
plane (along with some regions already excluded by some initial pioneering 
experiments) and the theoretical expectations of invisible axion models. The 
hope is that when both the Livermore and Kyoto experiments are completed, in 
3-5 years, one will know whether axions are, or are not, an important component 
of the dark matter in the Universe 


8 Perspectives on Dark Matter 


It is useful at this stage to try to bring some perspective on the issue of dark 
matter from a particle physics point of view. As we saw, particle physics pro¬ 
vides an interesting array of dark matter candidates. Among these, it appears 
that perhaps the neutralino LSP is the particle physics relic which is the most 

“Standard” in this context means invisible axions with an initial misallignment angle 
5 ~ 0(1) and ones where coherent axion oscillations dominate the energy density contribution. 
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plausible dark matter candidate. In the simplest supersymmetric extension of 
the standard model, the MSSM, there is a rather large range in parameter space 
which gives rise to a neutralino LSP that has ^ 0{1). In contrast, both for 
axions, gravitinos and neutrinos, the critical density in the Universe obtains only 
for some specific values of the parameters characterising these excitations (e.g. 
for axions one needs the scale of ?7(1)pq breaking, /, to be of 0(10^^ GeV)). 

Although, on the face of it, the above argument seems very reasonable, I 
am not sure it is totally compelling. For instance, in a similar vein one could 
argue also that having Ub = 0.05 is unnatural, since it requires a peculiar 
tuning of the nucleon mass! I believe a more sensible point of view to take is 
the following. Of all the cosmological scenarios, the inflationary scenario for the 
Universe appears to make the most sense. If this scenario is correct, then U = 1 
is a boundary condition one should seriously impose as a constraint on the sum 
of all the particle species which are important in the Universe today. That is, 
we should demand that^ 



(77) 


The particular weight of each of the components in Eq. (77) is a reflection 
of intrinsic particle physics properties. The only cosmological constraint is that 
the sum of the fli must add up to unity. So, if particle physics arguments lead to 
/ ~ 10^^ GeV, or ~ 5 eV, then that particular component will be important 
in the sum appearing in Eq. (77). From this point of view, “what you see is 
what you get”! If the parameters in the neutrino sector lead to some neutrino 
masses being in the eV range, then is an important component of U. If that 
is not the case, then is not important. So, from this view point, there is 
no difference in pedigree between dark matter which is a significant component 
for a range of particle physics parameters (like the LSP), or relics which are 
important only for the specific value of some particle physics parameters (like a 
KeV gravitino). 

Adopting this point of view then, it is perfectly sensible to have various 
particle physics excitations (say: baryons, neutralinos and neutrinos) play an 
important role in the Universe now. This is a welcome result, which is reinforced 
by the power spectrum of density fluctuations in the Universe. This spectrum 
also suggests that there is more than one component which contributes to the 
energy density of the Universe. Indeed, present data on this spectrum seems to 
be best fit by having a variety of matter components contributing. For instance, 
recent work by Primack and collaborators!^ suggests that the power spectrum 
of density fluctuations is optimally fit by having 


(78) 


Ub = 0.05; UcDM = 0.75; Uhdm = 0.20 . 


principle, one of the Oi could be the contribution from a cosmological constant. 
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These results are particularly interesting since the existence, or not, of HDM 
provides a critical constraint (arising from cosmology) on the particle physics 
which determines the neutrino mass matrix.^ If it were really possible to estab¬ 
lish the need for neutrino hot dark matter, for example through the influence it 
has on the angular spectrum of the CMBR,p^ then this, along with the con¬ 
straints imposed by neutrino oscillation experiments would do much to fix the 
shape of the neutrino mass spectrum. As we discussed earlier, if one can estab¬ 
lish both the need for neutrino hot dark matter (which necessitates probably 
that rrivi — (5 — 6 eV)) and of neutrino oscillations with small mass squared 
differences, then one is forced into a world of nearly degenerate neutrino masses, 
with rriv^ ~ 1 — 2 eV. Such a result would provide a compelling argument for 
renewing the direct searches in tritium beta decay for electron neutrino masses 
in the eV range. 


9 The Sakharov Conditions for Baryogenesis 


In a classic paper, in 1967, Andrei Sakharov[Q discussed the conditions neces¬ 
sary to obtain dynamically an asymmetry between matter and antimatter in the 
Universe. Sakharov’s conditions for obtaining this asymmetry are three-fold: |Q| 

(i) The underlying physical theory must possess processes that violate baryon 

number (B is not conserved). 

(ii) The interactions which lead to B-violation, in addition must violate C and 

CP. 

(iii) To establish this asymmetry dynamically, furthermore, the B-violating 
processes must be out of equilibrium in the Universe. 

Let me comment briefly on each of these points. First, it is pretty clear that if B 
is conserved then the total number of baryons minus anti-baryons is a constant 
in time. In this case, then the difference ne — ng is a constant that is set by 
some initial boundary conditions. Thus rj is not generated dynamically, but is 
just a reflection of these initial boundary conditions and one is left to wonder 
why one has a value rj ~ 

Similarly, it is also quite understandable why the second Sakharov condition 
is needed. If C and CP are good symmetries, one can transform ub into by 
one of these symmetry transformations. Hence, even if B were to be violated, 

^^This information, of course, is of relevance for experiments looking for neutrino 
oscillations. 
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but if C or CP were to be good symmetries, then one could never obtain a 
non-vanishing value for rj. 

The third Sakharov condition is slightly more subtle, but is also readily un¬ 
derstandable physically. Roughly speaking, B-violating decays serve to create 
a matter-antimatter asymmetry. However, this asymmetry is destroyed by in¬ 
verse decays. In thermal equilibrium, the rates for B-violating decays and their 
inverses are the same, hence ub — ug =0. 

It is useful to demonstrate this last fact explicitly. The rate of change of 
Aub = tt-b — ng as a result of B-violating processes, if these processes are in 
equilibrium, is given by the thermodynamic equation 

■ (79) 

Here 7 ^ is the rate of B-violation per unit volume and ^ is the chemical po¬ 
tential. At high temperatures, one can expand the exponential factors and the 
above expression reduces to 


c?Ang 2/i 

dt ~ ■ 


However, in this temperature regime, one has simply that 

4 o 

Ang = ■ 

TT^ 


(80) 


(81) 


Hence 


where Tjg 
(82) that 


dAuB 

dt 


TT 

'Y 


7R 

T3 


Ang = -- 


-Tj^AriB , 


(82) 


is just the rate for B-violation, since V = T Thus, it follows from 


Ang = (AnB)o exp 




(83) 


Eq. (83) tells one that, if B-violating processes are ever in equilibrium, then 
these processes serve to destroy any pre-existing asymmetry (AnB)o- This is a 
very nice result)^ since it tells us that the value of 77 one computes dynamically, 
as a result of B-violating processes going out of equilibrium, is independent of 
any initial asymmetry (AnB)o- Hence, the observed value of t] in the Universe 
now depends only on the B-violating (and C- and CP-violating) dynamics—due 
to particle physics—and on the cosmology which drives these processes out of 
equilibrium in the early Universe. 


I examine next cosmological circumstances (along with the relevant particle 
physics) which can lead to baryogenesis. 
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10 Baryogenesis at the GUT Scale: Issues and 
Challenges 


Grand Unified Theories (GUTs) were the first theories which explicitly real¬ 
ized Sakharov’s conditions for baryogenesis. 0 These theories naturally con¬ 
tain B-violating processes which also violate G and CP. An example is pro¬ 
vided by S't/(5),[^ in which the fermions of each generation are members of 
a 5 ^ (d£;eLt'L) and a 10 ^ (ui^di^; d^) representation^ and the ordinary 
Higgs doublet is augmented by a Higgs triplet y into a field 5h 

(x; —with X transforming under SU{S) x SU(2) xU{l) as X'^ (3, l)-i/ 3 . 

In SU{5), the Higgs quintet 5h can couple to the fermions in two separate ways 
[5h 5 10 and 5h 10 10], with the corresponding complex Yukawa couplings being 
sources for C and CP violation. These couplings allow the triplet Higgs field x 
to decay to both the dtx (B = 1/3) and ud (B = -2/3) final states. Hence, in 
SU (5) baryon number is clearly not conserved. 

Because one knows experimentally that baryon number is conserved to high 
accuracy,^ one knows that a theory like 5'C/(5), where the 5'C/(3), SU{2) and 
17(1) forces are unified, must break down to SU{3) x SU{2) x U{1) at a very 
high scale: Mx ~ 10^® — 10^® GeV.|^ This unification scale Mx is quite 
near the Planck scale Mp ~ 10^® GeV. We know that at temperatures near the 
Planck scale, T ~ Mp, the Universe is expanding very rapidly. Thus it is not 
surprising that the C, CP and B-violating decays of GUTs have rates which are 
slow with respect to the expansion rate of the Universe at T ~ Mx ■ That is 

^ ^ ^ re-viol. (T ~ Mx) ■ (84) 

Hence, the processes alluded above in GUTs also fulfill Sakharov’s third condition— 
that the relevant B-, C-, and CP-violating interactions be out of equilibrium in 
the early period of expansion of the Universe after the Big Bang. 

These qualitative features, however, in practice do not lead to successful 
simple scenarios for baryogenesis at the GUT scale. Although it is possible to 
obtain rj ^ 10“^° in some GUT models, these models have a number of generic 
difficulties which are worth discussing here. Again, it is useful to consider the 
SU{5) example alluded above to help focus on the source of these difficulties. 

In SU{5), the ratio ry is generated through the out of equilibrium decay of 

is convenient to describe all states in terms of how their left-handed components trans- 
form, using that '0£. 

^^The PDG[p^ gives a bound for the B-violating decay p 7r*^e+ of t{p —> 7r*^e"*") > 5 x 10^^ 
years. 
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the Higgs triplet x at temperatures T ^ Mx- One has 




( 85 ) 


Here H is a kinematical/dynamical factor related to the way the x decays go 
out of equilibrium, while is the baryon asymmetry proper: 


Ai?x = E 
/ 


Bf{T{x ^ /) - r(x - /)} 


ptotal 


( 86 ) 


reflecting the differences in the weighted ratio of x and x decays into particular 
final states / with different baryon number It should be clear from the 

form of Eq. ( 86 ) that AB;^ vanishes if C or CP is conserved, since then r(x ^ 

/) = r(x ^ /). 

It is easy to see that, for the example in question, one has simply 


AB^ = r — r 


where 


r(x ^ du) 

T = - 

ptotal 

Eq. (87) has three characteristics: 


r = 


r(x ^ dv) 

ptotal 
X 


(87) 

( 88 ) 


(i) It vanishes if there is no C or CP violation. This is obvious, since then 

f = r. 

(ii) ABjj. vanishes also if one includes only lowest order processes. Again this 
is easy to see since, at tree level, r = f. 

(iii) Finally, and less obviously, AB^ also vanishes if the underlying x-decays 
do not have an s-channel discontinuity. 


One can see these three conditions at work by examining schematically the 
contribution to r in the SU{5) model we discussed earlier. At one-loop order, 
these contributions are given by the graphs shown in Fig. 15. Let us denote by 70 
the rate associated with the tree graph decay in Fig. 15 and by 7 i/(M^ — ie) the 
contribution of the one-loop graph. In general both 70 and 71 are intrinsically 
complex as a result of the complex x couplings, while the dynamical quantity 
/(M^ — ie) has an imaginary part as a result of the associated loop integration. 

^®In the example discussed above f = dv or f = du. 


40 







10 



Figure 15: Graphs contributing to r at one-loop order 

A simple calculation, using Fig. 15, gives for the rate difference r — r the 
expression: 

(r-f) ~ |7o7 i/(M^- ie)|^ - |7o 

~ Im 7 o 7 i Im ie)) . (89) 

One sees that this rate difference vanishes unless there is both an intrinsic CP 
violating phase difference in the couplings involved [Im 7 o 7 i] as well as some 
imaginary part [Im (/(M^ ~*e)) ] arising from the (one-loop) scattering dynam¬ 
ics. In view of Eq. (89), one sees that the ratio rj = Ane/n-y is proportional 
to 

rj = AAB^ = A{r — f) ~ A Im 7o7i Im I . (90) 

The RHS of Eq. (90) embodies the essence of GUT baryogenesis. The 
ratio T] depends on the out of equilibrium dynamics [through A] and it van¬ 
ishes unless there is both an intrinsic CP and C violating phase [Im 707 *] 
and the GUT dynamics is rich enough to generate an s-channel discontinuity 
[Im I{M^ — ie)]. The knowledge of each of these individual pieces is clearly 
model-dependent and quite rudimentary, since we have no direct evidence for 
the existence of any GUTs ! Thus, at this stage, it is really not possible to 
deduce a firm prediction for rj. Even so, in general, one finds r] to be too small 
unless one further complicates the GUT dynamics. 

Let me illustrate the above point in the, by now familiar, 517(5) context. 
Without loss of generality one can make the 5 h 5 10 Higgs coupling matrix / 
real. Then it is easy to show that (for 3 families) the 5h 10 10 coupling matrix 
h has 3 phases. [Q So GUTs, because they involve further Higgs couplings, 
have more phases than the 3-family GKM phase connected with the couplings 
of the Higgs doublet $ to quarks. Even so, in this model, one cannot generate 
an intrinsic GP violating phase at one loop order. The tree and one-loop level 
contributions in Fig. 15, corresponding to the process 5 h ^ 10 5, give 

7o ~ / ; 7i ~ fhh^ ■ (91) 
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Figure 16: Interference graph giving r ^ 0 in SU{5) 


Hence 

Im 7o7i ~ Im Tr fhh^f^ = 0 . (92) 

One can check that other possible contributions to the decay 5 h ^ 10 5, 
involving gauge exchange in the t-channel rather than Higgs exchange, are sim¬ 
ilarly relatively real. As shown in Fig. 16, one can eventually |£l| obtain a 
non-vanishing ry for this model at higher order, from the interference of a tree- 
level process with a 3-loop process. The resulting however, 

AB^ ^ Im Tr[h^fphfph^h] (93) 

has such a large number of Yukawa couplings that f] is at best of 0(10“^®). @0 

This difficulty can be remedied by using more elaborate GUTs (or includ¬ 
ing more low-energy states). However, there are two further generic problems 
connected with these types of models which serve to dampen the enthusiasm 
for attributing baryogenesis in the Universe to some GUT processes. The first 
of these additional problems is related to monopoles. In general GUTs lead to 
an overproduction of monopoles in the early Universe, badly violating one of 
the main features we know about the Universe now-namely that the present 
Universe’s energy density is near the critical density p ~ Pc-lQ 

one invokes a fourth generation of quarks and leptons, it is possible to boost up r} 
to the desired 0(10“^'^) level even in this simple model. 
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’t Hooft and Polyakov showed that monopoles always form when a sym¬ 
metry group breaks down to a subgroup containing a 17(1) factor-like the stan¬ 
dard model group. Hence, if GUTs exist, one expects that in the very early 
Universe at T ~ Mx, during the GUT phase transition, magnetic monopoles 
are formed. Generically, these GUT monopoles are superheavy, having a mass 
of order Mm ^ ~ 10^"^ GeV. During the GUT phase transition, domains of 

the broken phase of the GUT group form which are of typical size ^ ~ ^ ~ . 

The superheavy GUT monopoles physically correspond to topological knots be¬ 
tween these domains and hence have a density um This density is 

comparable to the photon density at this stage of the Universe 

nM(Tc) ~ ~ Tc ~ . (94) 

However, such a large monopole density is extremely problematic, because of 
the large mass of the GUT monopoles .[|| Indeed, from (94) one deduces that 

IImI ~MMn^| (95) 

I now ' I now ^ ' 

completely in contradiction with what we know! 

Inflation provides a resolution of the monopole problem by inflating expo¬ 
nentially the size of the domains—essentially reducing the monopole density 
to one per observable Universe. However, to re-establish ry in such a scenario, 
one has to reheat the Universe after the inflationary period to temperatures 
Treheat ~ lO^"* —10^® GeV, which is difficult to achieve. HP! Thus, the monopole 
problem, even if it is resolved by inflation, argues against GUT baryogenesis. 

There is another argument which also provides ammunition against the idea 
that the baryon asymmetry in the Universe was produced at the GUT scale. 
As we will discuss shortly in more detail, it turns out that quantum effects in 
the electroweak interactions can lead to the violation of total fermion number 
(B-l-L—violation). 0 In the middle 1980’s Kuzmin, Rubakov and Shaposhnikov 
(KRS)[^ argued that these (B-l-L)-violating processes, which are extremely 
weak at T = 0, could become strong enough at temperatures near the elec¬ 
troweak phase transition, T ^ Mw, to go back into equilibrium in the Universe. 
The return of (B-l-L)-violating processes into equilibrium in the Universe at 
T <C Mx serves to erase any (B-fL)-asymmetry produced in the Universe at 
temperatures of the order of the GUT scale, T ^ Mx- Hence, only a (B-L) 
asymmetry produced by GUTs survives to low temperatures. 

This consideration kills, for example, the baryon number asymmetry one 
imagined was produced in the SU{5) example discussed earlier. It is easy to 
check that for y-decays 

= AL^ = r — f , (96) 

^®At such temperature the number density of monopoles produced after reheating is heavily 
suppressed by a Boltzmann factor. 
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so that 


AriB-L = AriB — AriL = 0 . (97) 

Thus, as a result of the KRS mechanism, in this case no baryon asymmetry 
survives at low temperatures, even if such an asymmetry were to be generated 
at the GUT scale by the out of equilibrium decays of the Higgs triplet x- Of 
course, one can invent more elaborate GUTs scenarios in which at the GUT 
scale one produces both a (B+L)- and a (B-L)-asymmetry, thereby bypassing 
this conundrum. 


11 The KRS Mechanism and Baryogenesis at 
the Electroweak Scale 


In this section I want to discuss further the KRS mechanism|^ because, be¬ 
sides erasing any previous (B-l-L)-asymmetry, it is possible that through this 
mechanism one can actually produce the observed baryon asymmetry in the 
Universe during the electroweak phase transition. This is an exciting possibil¬ 
ity, and one that has received considerable attention in recent years. In the 
Standard Model, both baryon number, B, and lepton number, L, are classical 
symmetries. That is, they are symmetries of the Standard Model Lagrangian:^® 

>CsM Asm • (98) 

However, because of the chiral nature of the electroweak interactions, at the 
quantum level both the baryon number current, Jg, and the lepton number 
current, Jg, are not conserved. Hence neither B, nor L, remains a good sym¬ 
metry at the quantum level, although their difference, B-L, is still a conserved 
quantum number. 

The violation of B and L in the standard model comes about as a result of the 
existence of chiral anomalies)^ in their respective currents. For our purposes, 
it suffices to focus only on the SU (2) gauge field contribution to this anomaly. 
The triangle graphs contributing to the anomalous divergence of Jg and Jg are 
shown in Fig. 17 and produce an equal divergence for both currents |8^ 

= -Ng^wrWa^. , (99) 

where Ng is the number of generations and a 2 = Glearly it follows then 

that 


neutrinos are massless, then the individual lepton numbers associated with electrons, 
muons and taus (Le, L^) are also SM Lagrangian symmetries. 
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qjarks 


leptons 


Figure 17: Triangle graphs contributing to the B and L anomalies 


( 100 ) 

These equations, per se, do not automatically lead to a violation of (B+L)- 
number. To get a change in B+L [A(B + L) ^ 0] requires having processes 
involving non Abelian gauge field configurations which have a non-trivial index 

V. 

^ ^ , ( 101 ) 

since, in view of (100), 

A(B + L) = 2Ag 1/. (102) 

’t HooftQ was the first to estimate the size of the amplitudes which contain 
gauge field configurations having such a non-trivial index u. These amplitudes 
arise in processes where the pure gauge field configurations at t = -l-oo and t = 
—oo differ by a so-called, “large” gauge transformation psf . In the A° = 0 gauge, 
pure gauge fields can be classified by how their associated gauge transformations 
go to unity at spatial infinity 


fl„(r) (103) 

One can show that the index v is related to the difference in the indices n± 
characterizing the gauge vacuum configurations at t = ±oo \v = n+ — n_].[Q 
’t Hooft’s estimate[|| of the size of the amplitudes leading to (B-l-L) violation 
essentially involved the WKB probability for tunneling from a vacuum charac¬ 
terized by index n to one where the index was n + u. His result Q| 

^(B-t-L)—violating ^ exp 

has a typical WKB form, involving the inverse of the gauge coupling con¬ 
stant squared in the exponent. However, since the weak coupling constant 


27r 

- u 

0.2 


(104) 
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squared a 2 is very small (02 ~ 1/30), the result (104) is extraordinarily tiny: 

^(B+L)-violating 10"®°" ! 


’t Hooft’s result (104) is valid at T = 0. What Kuzmin, Rubakov and 
ShaposhnikovQ realized was that the situation can be radically different in 
the early Universe, when the (B+L)-violating processes happen in a non-zero 
thermal background. When T ^ 0, the gauge vacuum change needed for 
A(B -|- L) ^ 0 transitions to happen can occur not only by tunneling, but also 
via a thermal fluctuation. In this latter case, the transition probability is not 
given by the square of the WKB amplitude (104), but instead by a Boltzman 
factor: 


violation ^ exp 


VojT) 

T 


(105) 


In the above, VoiT) is the (temperature dependent) height of the barrier which 
separates inequivalent gauge vacuum configurations. 


It turns out that one can estimate Vo also by semiclassical methods; in this 
case, by using a static solution of the electroweak theory with minimum energy 
and winding number n = 1/2. This solution, first found by Klinkhamer and 
Manton]^, has been dubbed by them a sphaleron. Essentially, one takes 
Vo{T) to be the energy associated with the sphaleron configuration in the pres¬ 
ence of a thermal bath: Vo{T) = ifsph(T). This energy has the typical form 
expected of a classical extended object. It is proportional to the mass of the 
gauge field associated with the symmetry which suffers the breakdown and is 
inversely proportional to the gauge coupling constant [c.f. the formula charac¬ 
terizing the monopole mass]. For the sphaleron, one has 

i?sph = ^^/(Mh/Mw) , (106) 

0:2 


where / is a function of order unity. 


Because the W mass, Mw{T), vanishes as the temperature T approaches 
the temperature of the electroweak phase transition, T —> Tew^ the probability 
of (B-l-L)-violating processes occurring in the Universe becomes large as the 
Universe’s temperature approaches Tew- This is basically the fundamental ob¬ 
servation made by Kuzmin, Rubakov and Shaposhnikov.|^ That is, one expects 
that 


-^(B-t-L)—violation(^) ~ CXp 


EsphjT) 

T 


"^EW 


(107) 


The original suggestion of KRS has been confirmed subsequently by much 
more detailed calculations. |Q Furthermore, one has found also a fast rate for 
(B-fL)-violation, above the temperature of the electroweak phase transition. [pl| 
These results are summarized below in a pair of formulas giving the transition 
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probability per unit volume, per unit time, for temperatures below and above 
the temperature of the electroweak phase transition. One finds: 


T(B+L) — violation 


T(B+L) — violation 


c 


Ml 


{a2Tf 


exp — 


Esp\i[T) 

T 


C'[a2Tf {T>Tew) 


{T < Tew) 


(108) 


where C and C" are constants of order one. These results imply that the rate 
of (B+L)-violating processes, originating in the standard model, is more rapid 
than the Universe’s expansion rate H ^ T^/Mp for rather a large temperature 
interval: 

r(B+L)-vioiation = > H for Tew ~ 10^ GeV < T < 1012 GeV . 

(109) 

A consequence of the above is that any (B+L)-asymmetry established above 
Tmax ^ 10i2 GeV (by, for example, some GUT processes) will get 

washed out. 


Given these results, two possibilities emerge for trying to explain the ob¬ 
served value of ry ~ 10“i°: 


i) The baryon-antibaryon asymmetry underlying ry is the result of a (B-L), or 

perhaps simply an L, asymmetry generated at high temperatures. Since 
the baryon number can be written as 

B=\{B + L) + ]^{B-L) , (110) 

and all the (B-l-L)-asymmetry is erased by the KRS mechanism, one needs 
to have some (B-L) asymmetry produced at high temperature to generate 
a non-vanishing value for 7y now. 

ii) The observed value of 7y is the result of processes occurring at the electroweak 

phase transition. Baryogenesis is simply the reflection of the violation of 
(B-fL) in the standard model-electroweak baryogenesis. 


For the remainder of this section, I want to discuss this latter possibility. 
This is a very intriguing suggestio n and one which has generated an enor¬ 
mous amount of interest recently. p6| It is clear that to be able to generate 7y 
at the electroweak phase transition, one needs this transition to be of first or¬ 
der, so as to get a deviation from thermal equilibrium. As one goes through 
the phase transition, the Higgs VEV, which vanished above Tew, jumps to a 
non-zero value {(piT*)) for temperatures below that of the electroweak phase 
transition. 
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Figure 18: Jump in the order parameter in units of the critical temperature 


However, to try to obtain through this non-equilibrium processes rj ^ 
needs much more than just having a first-order phase transition. There are 
actually two other main requirements. First, one must make sure that the 
asymmetry Atib-i-l created at the electroweak phase transition does not get 
erased by having the (B-l-L)-violating processes still be in equilibrium at T*. 
This requires that 


r(B-|-L)-violation(7"*) = C 


[ Kr 1 



exp 


Es^T*) 




< H{T*) 


rp*2 

Mp 


Numerically, this condition is equivalent to the requirement [86| 


EsphiT*)^,, WT*)) 

> 45 or ——- > 1 . 






( 111 ) 


( 112 ) 


That is, to avoid erasure of the produced Att-b+l after the electroweak phase 
transition, this transition must be strongly first order, giving rise to a large 
jump for the Higgs VEV. 


One can compute the jump in {(p(T*)) from the temperature-dependent ef¬ 
fective potential for the electroweak theory. Although there are a number of 
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uncertainties in this calculation,it is now generally agreed|Q that in the 
Standard Model with only one Higgs doublet one cannot get a sufficiently strong 
first-order phase transition, unless the Higgs boson is light. The present LEP 
bounds on the Higgs boson mass, Mh > 77.5 GeV|^ are already sufficiently 
strong so as to rule out this simplest version of the Standard Model as the 
source of the Universe’s matter-anti-matter asymmetry. This is made clear 
by Fig. 18, taken from|M , which shows that for Mh ~ 75 GeV one expects 


{(j)(T*))/T* ~ 0.5, in contradiction with the requirement of Eq. (112). For such 
high values of the Higgs mass, standard model processes could create a baryon 
asymmetry Atib+l at the electroweak phase transition, but this asymmetry 
would then get destroyed again at T*, since the rate of (B-l-L)-violation is still 
quite fast at this temperature compared to the Universe’s expansion rate. 


The situation, in this respect, is considerably better in supersymmetric ex¬ 
tensions of the Standard Model. Carena, Quiros and Wagner, Q for example, 
recently showed that the inequality (112) can be satisfied in the MSSM, pro¬ 
vided that one has a light stop, as well as small values for tan (3 and a reasonably 
light Higgs. The values for mi and mu required by |Q to allow for baryogenesis 
at the weak scale are low enough that these particles should be observable in 
the near future already at LEP 200 and/or at the Tevatron, when the Main In¬ 
jector is put into operation. So perhaps one may, in this way, soon get indirect 
evidence for electroweak baryogenesis. 

The second important requirement for believing that baryogenesis occurred 
at the electroweak scale is to actually be able to carry out a detailed dynamical 
calculation for 77, yielding 77 ~ 10 “^°. Even if the electroweak phase transition is 
sufficiently strongly first order, so that the produced Atib+l is not erased, it is 
not obvious that one can produce a big enough Atib+l so as to give 77 ^ 10“^°. 

In the case of electroweak baryogenesis, the matter-antimatter asymmetry 
is produced when bubbles of the true vacuum grow and fill up the Universe 
after the electroweak phase transition. As the Universe goes through the (as¬ 
sumed) first-order electroweak phase transition, CP-violating processes occur¬ 
ring in matter in the expanding bubble of true vacuum are crucial to estab¬ 
lishing a matter-antimatter asymmetry. The rate of (B-l-L)-violation per unit 
volume, 7(b-i-l)- violation) is rapid in the symmetric vacuum surrounding the bub¬ 
bles of true vacuum [7 °b+l)' 1 violation (ct2T)'^]. However, by assumption, if 
Eq. (112) holds, within the bubbles of true vacuum this same rate is negli¬ 
gible [7(6+^-violation — 0]- This rate difference is the key for establishing an 
asymmetry. 

A full dynamical calculation of 77 is very difficult. Nevertheless, one can 
sketch pictorially what is going on simply by thinking of the expanding true 
vacuum bubble as a wall sweeping through a plasma of quarks and antiquarks. 
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Because of CP-violating processes in the wall, the scattering of quarks and 
antiquarks off the wall, as this wall sweeps through the plasma, is not the 
same. This difference in scattering will create an excess, say, of antiquarks over 
quarks in the symmetric vacuum, —Auq, and an opposite excess of quarks over 
antiquarks, ArZg, in the true vacuum. The fast (B+L)-violating interactions 
in the symmetric vacuum, however, rapidly erase the antiquark excess, —Arig, 
leaving a quark excess, Arig, in the true vacuum bubble. This is the source of 
the baryon asymmetry. 


There are many issues one has to resolve, or understand, to perform a reli¬ 
able calculation of the above processes. For instance, what is the bubble wall 
thickness?; what is the bubble value velocity?; etc. Nevertheless, although an 
actual calculation is difficult, one can at least arrive at an order of magnitude 
estimate for 77 . Recall that the asymmetry is produced by the erasure of the 
antiquark excess in the symmetric vacuum. Hence 


V 


"TB+L violation 


T* 


T3 


sym vacuum 


C’a^T* . 


(113) 


If one assumes, quite naturally, that all time scales in the problem are set by 
1/T*,|Q then the proportionality constant in Eq. (113)—besides the factor of 
a 2 —is set by the amount of CP-violation in the quark and antiquark scattering 
off the bubble wall. Calling this factor ecp-vioiation, one arrives at the estimate 


77 — [ 0 ^ 2 ] CCP—violation — 10 CQP—violat 


(114) 


This estimate is confirmed by more detailed calculations, like those done recently 
by Huet and Nelson, 


Taking Eq. (114) as a reasonable guesstimate for electroweak baryogenesis, 
one sees that this process is effective in creating a sufficiently large baryon- 
antibaryon asymmetry only if there is enough CP-violation! That is, one needs 
to generate in the bubble walls at least ecp-vioiation 10“^. It turns out, 
however, that the standard model of flavor violation—the CKM model[p 8 ||— 
fails miserably in this task. In this case, because of the GIM mechanism, |]99{ 
there is no CP-violation unless the quark masses are different. Whence, one 
expects in this case that ecp-vioiation contains a number of GIM factors, which 
vanish if there is quark degeneracy. In particular, on general grounds, one 
expects 


[ecp- 


iCKM [A® sin ( 5 ] 2 2 \r 2 2 \/ 2 2^ 

violationj ~ “ "^c)(Wc “ mj 

■ {rnl - - ml){ml - ml)] . (115) 


The first square bracket contains the usual family mixing suppression factor 

Eq. (115) A is the sine of the Cabbibo angle and <5 is the CP-violating phase in the 
CKM model. 
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[ [100| while the second factor involves a product of GIM factors. The result 
which follows from Eq. (115), Ccp-vioiation ^10 is very small, falling far 
short of what is needed. 

Eq. (115), however, may be too naive an estimate. For example, one can 
avoid altogether the GIM suppression factor of Eq. (115) in models where 
there are some non-flavor violating sources of CP-violation at the electroweak 
scale. Examples of such models are provided by multi-Higgs models, or mod¬ 
els involving a supersymmetric extension of the Standard Model. Thus, if one 
believes that baryogenesis is really an electroweak-scale phenomenon, one is 
forced to contemplate theories which can provide a big enough ecp-vioiation- 
Particle physics theories which produce this result are necessarily enlargements 
of the Standard Model, since the Standard Model itself cannot provide enough 
ecp-vioiation- H baryogenesis occurred at the electroweak scale, these consid¬ 
erations argue that one is to expect both physics beyond the Standard Model, 
and other GP-violating phases besides the standard GKM phase, already at this 
scale. Gonversely, if one were to find these phenomena experimentally in the 
future, this would also provide indirect evidence for electroweak baryogenesis. 


12 Generating a B-asymmetry from an L-asymmetry 


If the matter-antimatter asymmetry is not generated by electroweak baryogene¬ 
sis then, as I mentioned earlier, this asymmetry must arise from processes which 
violate B-L at early times in the Universe. Because all (B-l-L)-asymmetries gen¬ 
erated above temperatures of order T ^ 10^^ GeV are erased by the KRS 
mechanism, purely L-violating processes effectively are equivalent to (B-L)- 
violating processes. However, these (B-L)-violating, or L-violating, processes 
cannot themselves go back into equilibrium after the asymmetry is generated, 
because that would serve again to erase this asymmetry. As we shall see, this 
last requirement has some (mild) implications for neutrino masses. 


I want to illustrate this last point by discussing briefly a specific model, due 
to Fukugita and Yanagida,[101 where what is violated is actually lepton num¬ 
ber. In the Fukugita-Yanagida scenario there are generic L-violating operators 
in the theory (arising from some GUT processes). These operators give rise 
both to a neutrino mass for i^lS and to lepton number violating processes. The 
simplest operator of this kind is one which involves the usual SU{2) x 17(1) 
left-handed lepton doublet L for the first generation and the Higgs field $: 


Cal=2 = -^L^CtL ■ + h.c. . 


(116) 


When $ is replaced by its VEV, this term gives rise to a mass term for the 
left-handed electron neutrino. At the same time, this term also contributes to 


51 



the L-violating process VeVe A straightforward calculation gives for the 

rate of L-violation 


FL-Violation = {na{VeVe ^ $$)) — ■ 


(117) 


If the masses of the neutrinos were to be large, this rate could be actually 
faster than the Universe’s expansion rate at temperatures below T ^ 10^^ GeV, 
where (B+L)-violating processes, due to the KRS mechanism|^, are themselves 
fast. If this were to be the case, then no matter-antimatter asymmetry 
would ever be generated at all! The necessary out of equilibrium condition 
bL-vioiation < H ^ T^/Mp, imposcs therefore a constraint on how larg e neutrino 
masses can be. Using the above result, one arrives at the bound] 1021 




< 


0.4 eV 


[T/1012 GeV] 


(118) 


That is, if neutrino masses are larger than this, then fast L-violating processes 
(in conjunction with the KRS mechanism) can erase any previously established 
matter-antimatter asymmetry. Of course, this bound is really rather soft in 
that it originates from only one possible type of L-violating interaction. Nev¬ 
ertheless, it is representative of a class of generic bounds which exist if one 
does not attribute the observed matter-antimatter asymmetry to electroweak 
processes. |l02| 


This said, however, one should mention that it is possible to avoid the above 
constraint. The simplest way to do this] 103] is to actually generate the matter- 
antimatter asymmetry of the Universe from L-violating processes which go out 
of equilibrium much below T ~ 10^^ GeV. This generally necessitates the 
introduction of right-handed neutrinos, with the out of equilibrium decays of pr 
generating the required asymmetry. Th e diffi culty in these scenarios, however, 
is producing a big enough asymmetry. [ 103 I will not discuss this matter in 
detail here. Suffice it to say that some successful models exist. These models 
have the peculiar feauture that rj is driven by the CP-violating phases in the 
neutrino sector! This last fact is easy to understand since these phases, 6^. are 
the ones which drive the lepton asymmetry [Anp ~ sin 6,^] and through the KRS 
mechanism, Anp = Anp. 


13 The Lessons of Baryogenesis for Particle Physics 


The above discussion of baryogenesis, either as a result of GUT models or 
through Standard Model processes has been very speculative. I believe, how¬ 
ever, that there are two overarching lessons one can draw from it. The first of 
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these is that there are a plethora of particle physics scenarios which can serve 
to generate a matter-antimatter asymmetry in the Universe. Thus, Sakharov’s 
intuition, that this asymmetry is dynamically generated and not the result of 
some peculiar initial boundary condition, is most likely true. 


The second lesson one draws from these disquisitions is that what is crucial 
for the whole issue of baryogenesis is the existence of other CP-violating phases, 
besides the usual CKM phase. An important goal, therefore, from a particle 
physics point of view is to try to discover these phases experimentally. This is 
a difficult, but perhaps not impossible task. I would like to end these lectures 
by making a few remarks on this point, particularly as it concerns electroweak 
baryogenesis. 


Electroweak baryogenesis suggests the presence of flavor diagonal CP- 
violating phases. These phasess arise quite naturally in multi-doublet Higgs 
models and in supersymmetric extension of the Standard Model. These phases, 
in general, do not contribute significantly to CP-violating quantities which are 
sensitive to the CKM phase—like the CP-asymmetries in B decays to CP self¬ 
conjugate states. 1 1^1 However, they can give rather large contribution to some 
CP-violating parameters which are small in the CKM model, like the electric 
dipole moment of the neutron. For example, supersymmetric extensions of the 
Standard Model give an electron dipole moment of the neutron which is of 


order] 105 


drfl — 


10 ^“sini^suSY 

^susY(GeV) 


-ecm 


(119) 


In the above 0 susy and Msusy are generic SUSY phases and masses. One sees 
that for Msusy ~ 100 GeV, the present limits for the electric dipole moment of 
the neutron, of order c?„ < 10“^® ecm,Q requires </>susy < 10“^. In fact, t here 
is no real explanation why the SUSY violating phases should be so small.] 1061 
So, it is obviously very important that one should push the experimental limit 
for dn beyond 10“^® ecm, as a supersymmetric signal may just be lurking around 
the corner! Unfortunately, this is a very difficult task in practice as experiments 
may have already reached their ultimate sensitivity limit. 


Fortunately, dn is not the only quantity which is sensitive to flavor diagonal, 
CP-violating phases. Another interesting measurable quantity, which perhaps 
is experimentally more accessible, is the transverse muon polarization in the 
decays K'^ —> This quantity measures the T-violating correlation: 


(P^) = 


[Sf, ■ {Pf, X p^) . 


( 120 ) 


Although (py) is not a purely CP-violating signal, the final state interactions 
in this decay which could also produce a transverse polarization are negligibly 
small [(py)FSi < 10“®[^|]. Thus a measurement of this quantity should test 
CP-violation. 
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What is interesting about (p^) is that this quantity vanishes in the Stan¬ 


dard Model. 1108] However, if there are CP-violating effective scalar interactions 
arising from physics beyond the standard model, one can get a significantly 
large transverse muon polarization. For instance, Grossman |109| finds that CP- 
violating phases in the Higgs sector (in multi-Higgs models) which satisfy the 
present bound on (i„, give a bound on ( p^) < 10“^. Remarkably, the present 
bounds on (p^) are precisely at this level [|110| 


(p^) = (-3.1 ± 5.3) X 10 


-3 


( 121 ) 


There is an experiment underway at KEK at the moment which hopes to push 
the error on (p^) to perhaps as low as (5(p^) ~ 5 x lO”"*.]!!!] Whether this 
can be achieved remains to be seen. However, in the near term, perhaps this is 
the best chance for finding some non-CKM sources of CP-violation. This is an 
exciting and important discovery window. If found, a non-vanishing value for 
(pi^) would have profound implications, not only for particle physics but also 
for cosmology. 
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